PhD Chapter 2
Ecological Status calculation and analyses
This series of files compile analyses done during Chapter 2.
All analyses have been done with R 4.0.4.
Click on the table of contents in the left margin to assess a specific analysis.
Click on a figure to zoom it
1. Indicator results
1.1. Scatterplots
Here are the plots presenting the value of each indicator at shallow and deep stations.
1.2. Ecological Quality Status
When relevant, we calculated an Ecological Quality Ratio (EQR) as established by the WFD and MSFD (which varies between 0 and 1). This ratio is calculated with the following equation:
\[ EQR = \frac{V_{ind} - R_{bad}}{R_{good} - R_{bad}} \]
- \(V_{ind}\) is the value of an indicator at a certain location
- \(R_{bad}\) is the reference value for a “bad” status
- \(R_{good}\) is the reference value for a “good” status
This ratio is then classed into Ecological Quality Status (EQS) categories, where reference values and limits for class transitions are specific to each indicator. Five classes are typically described:
- bad (red #FF0000)
- poor (orange #FFA500)
- moderate (yellow #EEEE00)
- good (green #228B22)
- high (blue #0000EE)
We calculated this ratio using different indicators, in order to compare their efficiency and relevance.
AMBI
We defined class thresholds using the methods from Borja et al. (2000) and Muxika et al. (2005).
M-AMBI
We defined class thresholds using the method from Muxika et al. (2007).
BENTIX
We defined class thresholds using the method from Simboura & Zenetos (2002).
BOPA
We defined class thresholds using the method from Dauvin & Ruellet (2007).
2. Non-metric Multidimensional Scaling
Species densities have been transformed with a (log+1) operation. Green points are stations within the bay, red points are within the archipelago.
Total density
Total biomass
W-Statistic
Specific richness
Shannon index
Margalef index
Simpson index
Pielou evenness
Taxonomic diversity
Functional richness
Functional evenness
Functional divergence
AMBI
M-AMBI
BENTIX
BOPA
3. Statistical relationships
3.1. Between indicators
In this section, we study the statistical relationships between indicators calculated above.
3.1.1. Covariation
Several types of models were considered to explore relationships: linear, quadratic, exponential and logarithmic. The model with the highest \(R^{2}\) is presented on each plot.
⚠️ Only linear models were implemented for now, as there are some bugs with the calculation of the others.
Total density
Total biomass
W-Statistic
Specific richness
Shannon index
Margalef index
Simpson index
Pielou evenness
Taxonomic diversity
Functional richness
Functional evenness
Functional divergence
AMBI
M-AMBI
BENTIX
BOPA
3.1.2. Correlation
Correlations have been calculated with Spearman’s rank coefficients.
| N | B | W | S | H | margalef | lambda | J | delta | FR | FE | FD | AMBI | M_AMBI | BENTIX | BOPA | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| N | 1 | 0.364 | -0.359 | 0.769 | 0.053 | 0.252 | -0.203 | -0.663 | -0.44 | 0.802 | -0.327 | 0.366 | 0.148 | 0.373 | 0.046 | 0.252 |
| B | 0.364 | 1 | 0.349 | 0.426 | 0.233 | 0.356 | 0.15 | -0.109 | 0.082 | 0.504 | 0.258 | 0.294 | -0.415 | 0.483 | 0.226 | 0.337 |
| W | -0.359 | 0.349 | 1 | 0.131 | 0.623 | 0.526 | 0.675 | 0.59 | 0.707 | -0.011 | 0.673 | 0.112 | -0.125 | 0.368 | -0.039 | 0.147 |
| S | 0.769 | 0.426 | 0.131 | 1 | 0.58 | 0.759 | 0.236 | -0.333 | -0.05 | 0.873 | 0.034 | 0.407 | 0.111 | 0.798 | 0.021 | 0.445 |
| H | 0.053 | 0.233 | 0.623 | 0.58 | 1 | 0.811 | 0.889 | 0.461 | 0.59 | 0.381 | 0.577 | 0.215 | -0.002 | 0.778 | -0.104 | 0.298 |
| margalef | 0.252 | 0.356 | 0.526 | 0.759 | 0.811 | 1 | 0.613 | 0.222 | 0.476 | 0.583 | 0.412 | 0.253 | -0.083 | 0.855 | 0.062 | 0.411 |
| lambda | -0.203 | 0.15 | 0.675 | 0.236 | 0.889 | 0.613 | 1 | 0.701 | 0.754 | 0.085 | 0.651 | 0.078 | -0.023 | 0.502 | -0.179 | 0.083 |
| J | -0.663 | -0.109 | 0.59 | -0.333 | 0.461 | 0.222 | 0.701 | 1 | 0.862 | -0.411 | 0.544 | -0.151 | -0.296 | 0.055 | 0.115 | -0.158 |
| delta | -0.44 | 0.082 | 0.707 | -0.05 | 0.59 | 0.476 | 0.754 | 0.862 | 1 | -0.195 | 0.514 | -0.109 | -0.279 | 0.263 | 0.065 | -0.099 |
| FR | 0.802 | 0.504 | -0.011 | 0.873 | 0.381 | 0.583 | 0.085 | -0.411 | -0.195 | 1 | -0.066 | 0.314 | 0.098 | 0.636 | 0.001 | 0.368 |
| FE | -0.327 | 0.258 | 0.673 | 0.034 | 0.577 | 0.412 | 0.651 | 0.544 | 0.514 | -0.066 | 1 | 0.343 | -0.223 | 0.395 | 0.01 | 0.207 |
| FD | 0.366 | 0.294 | 0.112 | 0.407 | 0.215 | 0.253 | 0.078 | -0.151 | -0.109 | 0.314 | 0.343 | 1 | -0.206 | 0.431 | 0.353 | 0.07 |
| AMBI | 0.148 | -0.415 | -0.125 | 0.111 | -0.002 | -0.083 | -0.023 | -0.296 | -0.279 | 0.098 | -0.223 | -0.206 | 1 | -0.318 | -0.781 | -0.184 |
| M_AMBI | 0.373 | 0.483 | 0.368 | 0.798 | 0.778 | 0.855 | 0.502 | 0.055 | 0.263 | 0.636 | 0.395 | 0.431 | -0.318 | 1 | 0.253 | 0.53 |
| BENTIX | 0.046 | 0.226 | -0.039 | 0.021 | -0.104 | 0.062 | -0.179 | 0.115 | 0.065 | 0.001 | 0.01 | 0.353 | -0.781 | 0.253 | 1 | 0.146 |
| BOPA | 0.252 | 0.337 | 0.147 | 0.445 | 0.298 | 0.411 | 0.083 | -0.158 | -0.099 | 0.368 | 0.207 | 0.07 | -0.184 | 0.53 | 0.146 | 1 |
| N | B | W | S | H | margalef | lambda | J | delta | FR | FE | FD | AMBI | M_AMBI | BENTIX | BOPA | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| N | 0 | 0.06769 | 0.07166 | 4.428e-06 | 0.7984 | 0.2149 | 0.3197 | 0.0002201 | 0.02459 | 8.389e-07 | 0.1027 | 0.06626 | 0.4694 | 0.06021 | 0.8253 | 0.2149 |
| B | 0.06769 | 8.884e-08 | 0.08109 | 0.0301 | 0.2512 | 0.07501 | 0.4626 | 0.597 | 0.6882 | 0.008604 | 0.2029 | 0.1453 | 0.03617 | 0.01241 | 0.2666 | 0.09228 |
| W | 0.07166 | 0.08109 | 8.884e-08 | 0.5234 | 0.0008811 | 0.006435 | 0.000222 | 0.001505 | 8.24e-05 | 0.959 | 0.0001671 | 0.5866 | 0.5398 | 0.06469 | 0.8513 | 0.4738 |
| S | 4.428e-06 | 0.0301 | 0.5234 | 2.315e-189 | 0.001917 | 7.044e-06 | 0.2459 | 0.09682 | 0.8101 | 6.131e-09 | 0.8675 | 0.0389 | 0.5902 | 1.037e-06 | 0.9202 | 0.02268 |
| H | 0.7984 | 0.2512 | 0.0008811 | 0.001917 | 8.884e-08 | 2.29e-06 | 1.665e-06 | 0.01789 | 0.001851 | 0.05487 | 0.002051 | 0.2906 | 0.9947 | 2.858e-06 | 0.6143 | 0.1388 |
| margalef | 0.2149 | 0.07501 | 0.006435 | 7.044e-06 | 2.29e-06 | 8.884e-08 | 0.001103 | 0.2759 | 0.01484 | 0.00176 | 0.03665 | 0.2123 | 0.6857 | 2.628e-08 | 0.7638 | 0.03696 |
| lambda | 0.3197 | 0.4626 | 0.000222 | 0.2459 | 1.665e-06 | 0.001103 | 8.884e-08 | 6.646e-05 | 1.509e-05 | 0.6792 | 0.000319 | 0.7063 | 0.9122 | 0.008977 | 0.3818 | 0.6856 |
| J | 0.0002201 | 0.597 | 0.001505 | 0.09682 | 0.01789 | 0.2759 | 6.646e-05 | 0 | 1.484e-08 | 0.03721 | 0.004045 | 0.4626 | 0.1419 | 0.7906 | 0.5771 | 0.44 |
| delta | 0.02459 | 0.6882 | 8.24e-05 | 0.8101 | 0.001851 | 0.01484 | 1.509e-05 | 1.484e-08 | 8.884e-08 | 0.34 | 0.007193 | 0.5947 | 0.1677 | 0.1943 | 0.7511 | 0.6316 |
| FR | 8.389e-07 | 0.008604 | 0.959 | 6.131e-09 | 0.05487 | 0.00176 | 0.6792 | 0.03721 | 0.34 | 0 | 0.7474 | 0.1183 | 0.6346 | 0.0004733 | 0.996 | 0.06463 |
| FE | 0.1027 | 0.2029 | 0.0001671 | 0.8675 | 0.002051 | 0.03665 | 0.000319 | 0.004045 | 0.007193 | 0.7474 | 0 | 0.08591 | 0.2736 | 0.0456 | 0.9616 | 0.3113 |
| FD | 0.06626 | 0.1453 | 0.5866 | 0.0389 | 0.2906 | 0.2123 | 0.7063 | 0.4626 | 0.5947 | 0.1183 | 0.08591 | 0 | 0.313 | 0.02783 | 0.0765 | 0.7352 |
| AMBI | 0.4694 | 0.03617 | 0.5398 | 0.5902 | 0.9947 | 0.6857 | 0.9122 | 0.1419 | 0.1677 | 0.6346 | 0.2736 | 0.313 | 8.884e-08 | 0.113 | 2.494e-06 | 0.3682 |
| M_AMBI | 0.06021 | 0.01241 | 0.06469 | 1.037e-06 | 2.858e-06 | 2.628e-08 | 0.008977 | 0.7906 | 0.1943 | 0.0004733 | 0.0456 | 0.02783 | 0.113 | 0 | 0.2133 | 0.005361 |
| BENTIX | 0.8253 | 0.2666 | 0.8513 | 0.9202 | 0.6143 | 0.7638 | 0.3818 | 0.5771 | 0.7511 | 0.996 | 0.9616 | 0.0765 | 2.494e-06 | 0.2133 | 0 | 0.4759 |
| BOPA | 0.2149 | 0.09228 | 0.4738 | 0.02268 | 0.1388 | 0.03696 | 0.6856 | 0.44 | 0.6316 | 0.06463 | 0.3113 | 0.7352 | 0.3682 | 0.005361 | 0.4759 | 0 |
| N | B | W | S | H | margalef | lambda | J | delta | FR | FE | FD | AMBI | M_AMBI | BENTIX | BOPA | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| N | 1 | -0.037 | -0.306 | 0.578 | -0.044 | 0.207 | -0.226 | -0.665 | -0.39 | 0.347 | -0.547 | 0.109 | 0.156 | 0.171 | -0.102 | 0.022 |
| B | -0.037 | 1 | 0.35 | 0.009 | -0.116 | -0.002 | -0.134 | -0.128 | -0.133 | 0.188 | 0.167 | -0.09 | -0.119 | -0.006 | 0.077 | 0.163 |
| W | -0.306 | 0.35 | 1 | 0.372 | 0.747 | 0.609 | 0.754 | 0.627 | 0.694 | 0.319 | 0.419 | -0.323 | -0.289 | 0.635 | -0.025 | -0.074 |
| S | 0.578 | 0.009 | 0.372 | 1 | 0.672 | 0.902 | 0.466 | -0.057 | 0.28 | 0.711 | -0.193 | -0.267 | -0.097 | 0.792 | -0.06 | -0.043 |
| H | -0.044 | -0.116 | 0.747 | 0.672 | 1 | 0.855 | 0.956 | 0.635 | 0.807 | 0.462 | 0.205 | -0.386 | -0.251 | 0.873 | -0.049 | -0.154 |
| margalef | 0.207 | -0.002 | 0.609 | 0.902 | 0.855 | 1 | 0.704 | 0.288 | 0.566 | 0.666 | 0.064 | -0.371 | -0.229 | 0.889 | -0.03 | -0.078 |
| lambda | -0.226 | -0.134 | 0.754 | 0.466 | 0.956 | 0.704 | 1 | 0.794 | 0.893 | 0.326 | 0.308 | -0.41 | -0.284 | 0.76 | 0.001 | -0.195 |
| J | -0.665 | -0.128 | 0.627 | -0.057 | 0.635 | 0.288 | 0.794 | 1 | 0.878 | -0.062 | 0.586 | -0.386 | -0.313 | 0.385 | 0.075 | -0.177 |
| delta | -0.39 | -0.133 | 0.694 | 0.28 | 0.807 | 0.566 | 0.893 | 0.878 | 1 | 0.21 | 0.432 | -0.498 | -0.297 | 0.603 | 0.019 | -0.224 |
| FR | 0.347 | 0.188 | 0.319 | 0.711 | 0.462 | 0.666 | 0.326 | -0.062 | 0.21 | 1 | -0.034 | -0.281 | -0.14 | 0.575 | -0.075 | 0.027 |
| FE | -0.547 | 0.167 | 0.419 | -0.193 | 0.205 | 0.064 | 0.308 | 0.586 | 0.432 | -0.034 | 1 | -0.118 | -0.199 | 0.093 | -0.097 | -0.134 |
| FD | 0.109 | -0.09 | -0.323 | -0.267 | -0.386 | -0.371 | -0.41 | -0.386 | -0.498 | -0.281 | -0.118 | 1 | 0.323 | -0.396 | -0.237 | 0.081 |
| AMBI | 0.156 | -0.119 | -0.289 | -0.097 | -0.251 | -0.229 | -0.284 | -0.313 | -0.297 | -0.14 | -0.199 | 0.323 | 1 | -0.52 | -0.697 | 0.115 |
| M_AMBI | 0.171 | -0.006 | 0.635 | 0.792 | 0.873 | 0.889 | 0.76 | 0.385 | 0.603 | 0.575 | 0.093 | -0.396 | -0.52 | 1 | 0.174 | -0.097 |
| BENTIX | -0.102 | 0.077 | -0.025 | -0.06 | -0.049 | -0.03 | 0.001 | 0.075 | 0.019 | -0.075 | -0.097 | -0.237 | -0.697 | 0.174 | 1 | -0.149 |
| BOPA | 0.022 | 0.163 | -0.074 | -0.043 | -0.154 | -0.078 | -0.195 | -0.177 | -0.224 | 0.027 | -0.134 | 0.081 | 0.115 | -0.097 | -0.149 | 1 |
| N | B | W | S | H | margalef | lambda | J | delta | FR | FE | FD | AMBI | M_AMBI | BENTIX | BOPA | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| N | 0 | 0.7422 | 0.005161 | 1.332e-08 | 0.6963 | 0.06213 | 0.04084 | 9.568e-12 | 0.0002929 | 0.001393 | 1.058e-07 | 0.3297 | 0.1605 | 0.1252 | 0.3631 | 0.8463 |
| B | 0.7422 | 0 | 0.001335 | 0.9369 | 0.2967 | 0.9868 | 0.2305 | 0.2524 | 0.2347 | 0.09033 | 0.1325 | 0.4203 | 0.2875 | 0.9597 | 0.4911 | 0.1432 |
| W | 0.005161 | 0.001335 | 0 | 0.0005736 | 0 | 1.266e-09 | 0 | 0 | 0 | 0.003636 | 0.000103 | 0.003214 | 0.008567 | 1.538e-10 | 0.8235 | 0.5089 |
| S | 1.332e-08 | 0.9369 | 0.0005736 | 0 | 4.887e-12 | 6.707e-31 | 1.024e-05 | 0.6104 | 0.01092 | 6.969e-14 | 0.08275 | 0.01521 | 0.3857 | 7.871e-19 | 0.5898 | 0.7001 |
| H | 0.6963 | 0.2967 | 0 | 4.887e-12 | 0 | 1.464e-24 | 0 | 0 | 0 | 1.565e-05 | 0.06474 | 0.0003808 | 0.02266 | 1.193e-26 | 0.661 | 0.167 |
| margalef | 0.06213 | 0.9868 | 1.266e-09 | 6.707e-31 | 1.464e-24 | 0 | 1.684e-13 | 0.008584 | 3.017e-08 | 8.921e-12 | 0.5687 | 0.0006022 | 0.03854 | 7.53e-29 | 0.7859 | 0.4853 |
| lambda | 0.04084 | 0.2305 | 0 | 1.024e-05 | 0 | 1.684e-13 | 0 | 0 | 0 | 0.002896 | 0.005111 | 0.0001512 | 0.009821 | 1.166e-16 | 0.9929 | 0.07979 |
| J | 9.568e-12 | 0.2524 | 0 | 0.6104 | 0 | 0.008584 | 0 | 0 | 0 | 0.5822 | 1.402e-08 | 0.0003808 | 0.004198 | 0.0003498 | 0.5017 | 0.1119 |
| delta | 0.0002929 | 0.2347 | 0 | 0.01092 | 0 | 3.017e-08 | 0 | 0 | 0 | 0.05853 | 5.912e-05 | 2.681e-06 | 0.006821 | 2.002e-09 | 0.8653 | 0.04351 |
| FR | 0.001393 | 0.09033 | 0.003636 | 6.969e-14 | 1.565e-05 | 8.921e-12 | 0.002896 | 0.5822 | 0.05853 | 0 | 0.7584 | 0.01075 | 0.209 | 1.572e-08 | 0.5053 | 0.8123 |
| FE | 1.058e-07 | 0.1325 | 0.000103 | 0.08275 | 0.06474 | 0.5687 | 0.005111 | 1.402e-08 | 5.912e-05 | 0.7584 | 0 | 0.2922 | 0.07322 | 0.4079 | 0.3841 | 0.2302 |
| FD | 0.3297 | 0.4203 | 0.003214 | 0.01521 | 0.0003808 | 0.0006022 | 0.0001512 | 0.0003808 | 2.681e-06 | 0.01075 | 0.2922 | 0 | 0.00305 | 0.0002311 | 0.03228 | 0.469 |
| AMBI | 0.1605 | 0.2875 | 0.008567 | 0.3857 | 0.02266 | 0.03854 | 0.009821 | 0.004198 | 0.006821 | 0.209 | 0.07322 | 0.00305 | 0 | 5.386e-07 | 3.684e-13 | 0.3041 |
| M_AMBI | 0.1252 | 0.9597 | 1.538e-10 | 7.871e-19 | 1.193e-26 | 7.53e-29 | 1.166e-16 | 0.0003498 | 2.002e-09 | 1.572e-08 | 0.4079 | 0.0002311 | 5.386e-07 | 0 | 0.1174 | 0.3835 |
| BENTIX | 0.3631 | 0.4911 | 0.8235 | 0.5898 | 0.661 | 0.7859 | 0.9929 | 0.5017 | 0.8653 | 0.5053 | 0.3841 | 0.03228 | 3.684e-13 | 0.1174 | 0 | 0.1817 |
| BOPA | 0.8463 | 0.1432 | 0.5089 | 0.7001 | 0.167 | 0.4853 | 0.07979 | 0.1119 | 0.04351 | 0.8123 | 0.2302 | 0.469 | 0.3041 | 0.3835 | 0.1817 | 0 |
3.2. Between indicators and abiotic parameters
In this section, we study the statistical relationships between indicators calculated above and different abiotic parameters, in order to understand how well they can be used to detect perturbations.
3.2.1. Covariation
Several types of models were considered to explore relationships: linear, quadratic, exponential and logarithmic. The model with the highest \(R^{2}\) is presented on each plot.
⚠️ Only linear models were implemented for now, as there are some bugs with the calculation of the others.
Total density
Total biomass
W-Statistic
Specific richness
Shannon index
Margalef index
Simpson index
Pielou evenness
Taxonomic diversity
Functional richness
Functional evenness
Functional divergence
AMBI
M-AMBI
BENTIX
BOPA
3.2.2. Correlation
Correlations have been calculated with Spearman’s rank coefficients.
| N | B | W | S | H | margalef | lambda | J | delta | FR | FE | FD | AMBI | M_AMBI | BENTIX | BOPA | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| depth | -0.183 | 0.151 | 0.361 | 0.116 | 0.345 | 0.252 | 0.206 | 0.375 | 0.361 | -0.112 | 0.21 | 0.196 | -0.324 | 0.325 | 0.47 | 0.022 |
| om | -0.073 | 0.089 | 0.089 | 0.107 | 0.089 | 0.187 | -0.098 | 0.013 | 0.022 | -0.036 | -0.004 | 0.092 | -0.432 | 0.311 | 0.449 | 0.347 |
| gravel | 0.096 | 0.055 | -0.026 | -0.024 | 0.071 | -0.005 | 0.196 | 0.143 | 0.155 | 0.223 | 0.105 | -0.194 | 0.319 | -0.151 | -0.365 | -0.052 |
| sand | 0.245 | -0.092 | -0.016 | 0.229 | 0.104 | 0.082 | 0.095 | -0.118 | 0.009 | 0.234 | -0.108 | 0.252 | 0.471 | -0.051 | -0.332 | -0.273 |
| silt | -0.181 | 0.141 | 0.041 | -0.136 | -0.071 | -0.03 | -0.117 | 0.049 | -0.027 | -0.182 | 0.105 | -0.136 | -0.474 | 0.132 | 0.352 | 0.252 |
| clay | -0.179 | 0.045 | 0.021 | -0.202 | 0.048 | -0.181 | 0.168 | 0.213 | 0.046 | -0.066 | 0.085 | -0.26 | -0.119 | -0.122 | 0.008 | -0.023 |
| arsenic | -0.456 | -0.423 | -0.1 | -0.474 | -0.113 | -0.182 | 0.038 | 0.283 | 0.225 | -0.428 | -0.023 | -0.595 | -0.014 | -0.336 | -0.124 | 0.015 |
| cadmium | -0.076 | -0.416 | -0.401 | -0.196 | -0.182 | -0.152 | -0.174 | -0.066 | -0.059 | -0.146 | -0.301 | -0.156 | -0.171 | -0.11 | -0.013 | -0.077 |
| chromium | -0.522 | -0.591 | 0.022 | -0.359 | 0.033 | -0.053 | 0.06 | 0.287 | 0.225 | -0.498 | 0.031 | -0.351 | 0.006 | -0.141 | -0.019 | -0.017 |
| copper | -0.546 | -0.506 | 0.138 | -0.327 | 0.05 | -0.018 | 0.068 | 0.315 | 0.23 | -0.426 | 0.059 | -0.372 | 0.092 | -0.174 | 0.006 | 0.011 |
| iron | -0.49 | -0.391 | 0.145 | -0.309 | 0.054 | -0.006 | 0.075 | 0.269 | 0.262 | -0.474 | 0.083 | -0.287 | -0.061 | -0.103 | 0.063 | -0.037 |
| manganese | -0.306 | -0.531 | -0.109 | -0.276 | -0.047 | -0.073 | -0.017 | 0.092 | 0.099 | -0.459 | -0.037 | -0.269 | 0.016 | -0.12 | 0.028 | 0.041 |
| mercury | -0.517 | -0.355 | 0.063 | -0.392 | 0.039 | -0.113 | 0.109 | 0.423 | 0.26 | -0.499 | 0.132 | -0.379 | -0.112 | -0.119 | 0.056 | 0.021 |
| lead | -0.549 | -0.499 | 0.067 | -0.381 | 0.022 | -0.034 | 0.086 | 0.365 | 0.237 | -0.42 | 0.062 | -0.4 | -0.013 | -0.174 | 0.037 | 0.148 |
| zinc | -0.515 | -0.613 | 0.029 | -0.343 | 0.058 | -0.037 | 0.097 | 0.32 | 0.259 | -0.465 | -0.011 | -0.411 | 0.063 | -0.172 | -0.031 | -0.032 |
| N | B | W | S | H | margalef | lambda | J | delta | FR | FE | FD | AMBI | M_AMBI | BENTIX | BOPA | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| depth | 0.3718 | 0.461 | 0.06988 | 0.571 | 0.08426 | 0.2148 | 0.3138 | 0.05915 | 0.06988 | 0.5853 | 0.3039 | 0.3373 | 0.1065 | 0.1047 | 0.01535 | 0.916 |
| om | 0.7237 | 0.6661 | 0.6636 | 0.6042 | 0.6636 | 0.3587 | 0.6323 | 0.9484 | 0.9175 | 0.8605 | 0.9841 | 0.6537 | 0.02864 | 0.1222 | 0.02133 | 0.08199 |
| gravel | 0.6414 | 0.791 | 0.8982 | 0.9083 | 0.7293 | 0.9803 | 0.3377 | 0.4863 | 0.4493 | 0.2741 | 0.6109 | 0.3434 | 0.1122 | 0.461 | 0.06668 | 0.7995 |
| sand | 0.228 | 0.6549 | 0.9392 | 0.2613 | 0.6145 | 0.6902 | 0.6429 | 0.5648 | 0.9643 | 0.249 | 0.5987 | 0.2145 | 0.01519 | 0.8035 | 0.09799 | 0.1778 |
| silt | 0.3755 | 0.4914 | 0.8422 | 0.5064 | 0.7286 | 0.884 | 0.5683 | 0.8112 | 0.8958 | 0.3732 | 0.6092 | 0.5079 | 0.01437 | 0.5192 | 0.07773 | 0.2134 |
| clay | 0.3828 | 0.8286 | 0.918 | 0.3216 | 0.8163 | 0.3756 | 0.413 | 0.2965 | 0.8224 | 0.7472 | 0.68 | 0.1992 | 0.5637 | 0.5511 | 0.9705 | 0.9109 |
| arsenic | 0.01931 | 0.03116 | 0.6272 | 0.01449 | 0.5839 | 0.3727 | 0.8547 | 0.1608 | 0.2693 | 0.02914 | 0.9105 | 0.001346 | 0.9443 | 0.09292 | 0.5461 | 0.9406 |
| cadmium | 0.7107 | 0.03431 | 0.04234 | 0.337 | 0.3728 | 0.4574 | 0.3962 | 0.7489 | 0.7743 | 0.4763 | 0.1345 | 0.4471 | 0.4042 | 0.5918 | 0.9498 | 0.7077 |
| chromium | 0.006246 | 0.001481 | 0.9142 | 0.07201 | 0.8722 | 0.7971 | 0.7703 | 0.1547 | 0.2691 | 0.009577 | 0.88 | 0.07851 | 0.9762 | 0.4924 | 0.926 | 0.9362 |
| copper | 0.003903 | 0.00909 | 0.4983 | 0.1029 | 0.8098 | 0.9308 | 0.7407 | 0.1171 | 0.257 | 0.02999 | 0.7728 | 0.06128 | 0.6539 | 0.3961 | 0.9775 | 0.958 |
| iron | 0.01112 | 0.0489 | 0.4771 | 0.1247 | 0.7917 | 0.9787 | 0.7155 | 0.1831 | 0.1949 | 0.01445 | 0.6878 | 0.1548 | 0.7686 | 0.618 | 0.7587 | 0.859 |
| manganese | 0.1279 | 0.005904 | 0.5945 | 0.172 | 0.8202 | 0.723 | 0.9334 | 0.6549 | 0.6275 | 0.01828 | 0.8592 | 0.1837 | 0.9388 | 0.5603 | 0.8905 | 0.8406 |
| mercury | 0.006883 | 0.0762 | 0.7584 | 0.04779 | 0.8516 | 0.5806 | 0.5945 | 0.03148 | 0.1986 | 0.009534 | 0.5204 | 0.0563 | 0.5852 | 0.5626 | 0.7842 | 0.9176 |
| lead | 0.003659 | 0.009392 | 0.7461 | 0.05509 | 0.9141 | 0.8709 | 0.6754 | 0.06698 | 0.2429 | 0.03246 | 0.7638 | 0.04293 | 0.9484 | 0.3944 | 0.8577 | 0.4699 |
| zinc | 0.007047 | 0.0008679 | 0.8866 | 0.08579 | 0.7792 | 0.8566 | 0.6358 | 0.1105 | 0.2023 | 0.01656 | 0.959 | 0.03695 | 0.7601 | 0.4017 | 0.8813 | 0.8774 |
| N | B | W | S | H | margalef | lambda | J | delta | FR | FE | FD | AMBI | M_AMBI | BENTIX | BOPA | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| depth | -0.17 | -0.074 | 0.3 | 0.206 | 0.447 | 0.309 | 0.467 | 0.469 | 0.534 | 0.259 | 0.267 | -0.375 | -0.166 | 0.401 | 0.032 | -0.231 |
| om | -0.101 | -0.037 | -0.155 | -0.248 | -0.13 | -0.263 | -0.047 | 0.018 | -0.149 | -0.206 | -0.205 | 0.009 | -0.092 | -0.157 | 0.266 | 0.203 |
| gravel | -0.013 | 0.205 | 0.239 | 0.123 | 0.076 | 0.146 | 0.061 | 0.058 | 0.116 | 0.25 | 0.095 | -0.115 | -0.085 | 0.118 | -0.055 | -0.023 |
| sand | 0.004 | -0.015 | 0.194 | 0.126 | 0.179 | 0.181 | 0.142 | 0.111 | 0.232 | 0.045 | 0.178 | 0.034 | 0.06 | 0.15 | -0.263 | -0.309 |
| silt | 0.072 | -0.019 | -0.219 | -0.106 | -0.202 | -0.187 | -0.193 | -0.195 | -0.286 | -0.09 | -0.208 | 0.034 | -0.035 | -0.139 | 0.231 | 0.341 |
| clay | -0.019 | -0.015 | -0.021 | -0.004 | -0.007 | 0.011 | -0.013 | -0.06 | -0.044 | 0.126 | -0.092 | 0.011 | -0.034 | -0.024 | 0.041 | 0.104 |
| arsenic | -0.031 | -0.079 | -0.288 | -0.27 | -0.285 | -0.317 | -0.223 | -0.169 | -0.291 | -0.254 | -0.204 | 0.207 | -0.048 | -0.24 | 0.225 | 0.331 |
| cadmium | -0.079 | -0.077 | -0.294 | -0.317 | -0.289 | -0.333 | -0.233 | -0.154 | -0.322 | -0.317 | -0.184 | 0.283 | 0.04 | -0.295 | 0.137 | 0.276 |
| chromium | -0.078 | 0.034 | -0.267 | -0.311 | -0.329 | -0.357 | -0.266 | -0.174 | -0.353 | -0.287 | -0.221 | 0.198 | -0.026 | -0.295 | 0.235 | 0.36 |
| copper | -0.069 | 0.006 | -0.235 | -0.318 | -0.285 | -0.365 | -0.222 | -0.169 | -0.348 | -0.322 | -0.252 | 0.293 | -0.06 | -0.266 | 0.254 | 0.31 |
| iron | -0.228 | 0.072 | -0.266 | -0.449 | -0.363 | -0.454 | -0.28 | -0.094 | -0.353 | -0.363 | -0.189 | 0.182 | 0.035 | -0.379 | 0.098 | 0.328 |
| manganese | -0.079 | 0.049 | -0.258 | -0.308 | -0.313 | -0.36 | -0.254 | -0.163 | -0.344 | -0.285 | -0.275 | 0.192 | -0.03 | -0.277 | 0.196 | 0.376 |
| mercury | 0.019 | 0.05 | -0.203 | -0.217 | -0.251 | -0.278 | -0.22 | -0.212 | -0.329 | -0.283 | -0.267 | 0.287 | -0.027 | -0.205 | 0.231 | 0.297 |
| lead | -0.041 | -0.024 | -0.278 | -0.301 | -0.306 | -0.347 | -0.256 | -0.195 | -0.343 | -0.273 | -0.238 | 0.283 | 0.01 | -0.282 | 0.193 | 0.334 |
| zinc | -0.069 | -0.053 | -0.286 | -0.342 | -0.308 | -0.379 | -0.241 | -0.182 | -0.349 | -0.326 | -0.218 | 0.335 | 0.004 | -0.306 | 0.18 | 0.303 |
| N | B | W | S | H | margalef | lambda | J | delta | FR | FE | FD | AMBI | M_AMBI | BENTIX | BOPA | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| depth | 0.1261 | 0.511 | 0.006261 | 0.06316 | 2.597e-05 | 0.004791 | 9.909e-06 | 8.77e-06 | 2.435e-07 | 0.01902 | 0.0154 | 0.0005261 | 0.1352 | 0.0001889 | 0.7763 | 0.03669 |
| om | 0.3675 | 0.7435 | 0.1635 | 0.02459 | 0.2456 | 0.01699 | 0.6754 | 0.8707 | 0.1826 | 0.0639 | 0.06448 | 0.9389 | 0.4133 | 0.158 | 0.01573 | 0.06805 |
| gravel | 0.9053 | 0.06425 | 0.03055 | 0.2716 | 0.4958 | 0.1918 | 0.5857 | 0.6039 | 0.2993 | 0.02327 | 0.3945 | 0.3057 | 0.4482 | 0.2916 | 0.6212 | 0.8381 |
| sand | 0.9702 | 0.8968 | 0.08083 | 0.2612 | 0.1078 | 0.1033 | 0.2026 | 0.3188 | 0.03598 | 0.69 | 0.1086 | 0.7624 | 0.5932 | 0.1775 | 0.01718 | 0.004724 |
| silt | 0.5178 | 0.8642 | 0.04799 | 0.3426 | 0.06857 | 0.09176 | 0.08247 | 0.07956 | 0.009136 | 0.4216 | 0.06117 | 0.7603 | 0.7533 | 0.2128 | 0.03641 | 0.001735 |
| clay | 0.868 | 0.8947 | 0.8486 | 0.9691 | 0.9523 | 0.9206 | 0.9067 | 0.5947 | 0.6935 | 0.2609 | 0.411 | 0.922 | 0.7611 | 0.833 | 0.714 | 0.3521 |
| arsenic | 0.7815 | 0.4799 | 0.008696 | 0.01399 | 0.009455 | 0.003656 | 0.04401 | 0.1283 | 0.007928 | 0.02149 | 0.06595 | 0.06212 | 0.6656 | 0.02989 | 0.04228 | 0.002378 |
| cadmium | 0.4801 | 0.4899 | 0.00741 | 0.003699 | 0.008482 | 0.002217 | 0.03551 | 0.1676 | 0.003198 | 0.003684 | 0.09801 | 0.01005 | 0.7245 | 0.007097 | 0.2204 | 0.01204 |
| chromium | 0.4853 | 0.7593 | 0.01533 | 0.004395 | 0.002539 | 0.0009817 | 0.01583 | 0.1185 | 0.00114 | 0.008854 | 0.04601 | 0.07398 | 0.8178 | 0.007079 | 0.03379 | 0.0008835 |
| copper | 0.5384 | 0.9568 | 0.03393 | 0.00356 | 0.009581 | 0.0007434 | 0.04513 | 0.1294 | 0.001377 | 0.003164 | 0.02234 | 0.007541 | 0.5946 | 0.01584 | 0.02136 | 0.004539 |
| iron | 0.0391 | 0.5209 | 0.01611 | 2.357e-05 | 0.0008703 | 1.821e-05 | 0.01105 | 0.3981 | 0.001222 | 0.0008618 | 0.08931 | 0.1023 | 0.7543 | 0.0004477 | 0.3805 | 0.002633 |
| manganese | 0.4802 | 0.6632 | 0.01965 | 0.004816 | 0.004388 | 0.0008973 | 0.0214 | 0.1436 | 0.001634 | 0.009644 | 0.01279 | 0.08341 | 0.7877 | 0.01168 | 0.07685 | 0.0004928 |
| mercury | 0.8681 | 0.6525 | 0.06723 | 0.0502 | 0.02294 | 0.01157 | 0.0472 | 0.05604 | 0.002681 | 0.01032 | 0.01564 | 0.009027 | 0.8123 | 0.06473 | 0.03687 | 0.006739 |
| lead | 0.7116 | 0.8275 | 0.01159 | 0.005972 | 0.005188 | 0.001414 | 0.0201 | 0.07908 | 0.001615 | 0.01299 | 0.03097 | 0.01008 | 0.9308 | 0.01013 | 0.08155 | 0.002162 |
| zinc | 0.5373 | 0.6367 | 0.009233 | 0.001654 | 0.00481 | 0.0004453 | 0.02904 | 0.1017 | 0.001332 | 0.002824 | 0.04918 | 0.002117 | 0.9743 | 0.005254 | 0.1064 | 0.005722 |
4. Indicator robustness
In this section, we are calculating values of the indicators for shallow and deep stations using a bootstrap method, so that we have an idea of the robustness of each measure.
| True mean | Bootstrap Mean bias | Boostrap 95% CI | |
|---|---|---|---|
| S | 9.192 | -0.03569 | [9.1708;9.2852] |
| N | 138.7 | -0.3769 | [136.8829;141.2555] |
| B | 7.352 | -0.05717 | [7.1399;7.6789] |
| W | 0.0109 | -0.0164 | [0.0269;0.0277] |
| H | 1.353 | -0.006664 | [1.3537;1.3663] |
| margalef | 1.926 | -0.01379 | [1.9302;1.949] |
| lambda | 0.6202 | -0.002792 | [0.6205;0.6255] |
| J | 0.6564 | -0.00366 | [0.6572;0.663] |
| delta | 51.66 | -0.3568 | [51.7936;52.2472] |
| FR | 23.35 | -3.171 | [26.111;26.9268] |
| FE | 0.5542 | 0.002318 | [0.5495;0.5543] |
| FD | 0.7656 | -0.007396 | [0.7701;0.7759] |
| True mean | Bootstrap Mean bias | Boostrap 95% CI | |
|---|---|---|---|
| S | 13.99 | -0.009295 | [13.964;14.0302] |
| N | 89.13 | -0.1556 | [88.8068;89.7726] |
| B | 8.72 | 0.0243 | [8.5538;8.838] |
| W | 0.02522 | -0.007476 | [0.0325;0.0329] |
| H | 1.952 | 0.0002019 | [1.9489;1.9551] |
| margalef | 3.046 | -7.372e-05 | [3.0402;3.052] |
| lambda | 0.7699 | 0.0001827 | [0.7687;0.7707] |
| J | 0.7601 | 0.0003165 | [0.7588;0.7608] |
| delta | 63.48 | 0.01433 | [63.386;63.5518] |
| FR | 31.76 | -7.59 | [38.827;39.8772] |
| FE | 0.6324 | -0.001514 | [0.6331;0.6347] |
| FD | 0.8282 | 0.01091 | [0.8165;0.8181] |